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Abstract 

We present a method to analyze the metal-insulator transition (MIT) due to the band overlap 
mechanism. It is based on a model with the knowledge of the homogeneous electron gas, combined 
with results based on the quasiparticle self-consistent GW method. Because of the long-range 
nature of the Coulomb interaction, the MIT occurs as the first-order phase transition, that is, the 
band gap becomes negative (band overlap) suddenly at some critical lattice constant. 
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FIG. 1: (color online) Illustration of the first-order metal-insulator transition by the band overlap 
mechanism. The band gap Eg is gradually being reduced when the lattice constant a is getting 
smaller. Then it suddenly jumps to the metallic phase; the band overlap -EqvI appears as shown 
by the broken line. In our method explained in Figl2] around (see text), we assume the rigid band 
shifts. 



The metal-insulator transitions (MIT) is important not only for fundamental physics, 
but also for the high-pressure physics, or for its potential applicability to electrical or op- 
tical switches, e.g, YH3 Among the possible mechanisms of the MIT {2!, we focus on 
the band-overlap mechanism which is the simplest in the sense that it is described within 
the one-particle picture. Here we present a theoretical treatment at zero temperature with- 
out phonons. The MIT is illustrated in Fig{T] (we explain it below). Though the density 
functional theory (DFT) can give an one-particle picture represented by its Kohn-Sham 
eigenvalues and eigenfunctions, the eigenvalues can not be identified as the quasiparticle 
energies (QPEs); it is well known that the DFT predicts too small Eg. The problem is not 
in the local density approximation (LDA) usually used in DFT; as shown by the calculations 
with the optimized effective potential method with the exact exchange plus correlation in the 
random phase approximation (EXX+RPA), Kotani showed that the true Kohn-Sham eigen- 
values without LDA are only slightly larger than the eigenvalues in LDA for semiconductors 
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31]. It is confirmed by other groups recently 
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Thus it is necessary to use a method beyond DFT to obtain the reasonable quasiparticle 
(QP) picture, or rather the one-particle static Hamiltonian which represent the QP. For 
example, for MIT for bcc Hydrogen, GW calculations is used by Kioupakis, Zhang,Cohen, 
and Louie [s]. The GW approximation (GWA) [6] is theoretically reasonable to obtain 
QPEs, however, the reliability of the usual GWA is limited because it is the one-shot (per- 
turbative) calculations starting from the solution in LDA. This can cause a problem in cases; 
for example, we may get insulator solution in GWA starting from the metallic solution in 
LDA. In such a case, its reliability is questionable. This situation can be corrected by the self- 
consistent perturbation idea. In fact, they used a simplified version of the self-consistency 
to determine U in the LDA+U+GW calculation as was done by Aryasetiawan and Gun- 
narsson It is also used recently for the MIT problem for VO2 by Sakuma Miyake and 
Aryasetiawan [s^. However, such a self-consistency is only for a degree of freedom; a satisfac- 
tory version of self-consistency is formulated by the quasiparticle self-consistent GW method 
(QSGPV). Then all the degree of non-local static potential is determined self-consistently. 
We had shown that QSGVF reproduce the band gap Eg for wide range of semiconductors 
and insulators very well [q, [io|. Thus it is reasonable to apply QSGWto analyze the MIT 
for the band-overlap mechanism. 

The QS GW is taken as an approximation to the rigorous theory. To see it, note that 
the self-consistent procedure in QSGW jo], Q, 11 1 can be divided into two parts. One is the 
determination of the one-particle static Hamiltonian from the given Hartree potential 
plus self-energy Vh + S(r, r', u), the other is the perturbative calculation of Vh + S(r, r', u) 
in GWA starting from H^. The former is a recipe to extract which represent the QP; 
note that iterative procedure is required to determine even for fixed Vh + S(r,r',ct;) 
(see the norm minimization formalism P]). The latter can be replaced by the rigorous 
procedure at least as a thought experiment. If we were able to do it, we would have reached 
to the exact ground state with reasonable H'^. A practical method to improve QSGH^is 
presented by Shishkin, Marsman and Kresse [3], where the excitonic effects (correlational 
motion of the electron-hole) are included in the polarization function when we evaluate the 
screened Coulomb interaction W. They succeeded to give systematic improvement for the 
overestimation of Eg in the QSGW. 

Here we present a method to analyze the MIT based on the knowledge of homogeneous 
electron gas in combination with the insulator solutions given by QSGW (or its extension). 
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As seen in Fig{Tl we have valence and conduction bands with Eg for the insulator phase. Eg 
is monotonically getting smaller when the lattice constant a is being reduced by external 
pressure. If we have a first-order MIT, Eg suddenly jumps from some positive value to the 
negative value (band overlap) at the critical lattice constant a = a^. This occurs if the 
total energy for the metal phase becomes lower than that of the insulator. Then we see 
electron pocket in valence band and hole pocket in conduction band, which are specified by 
the Fermi energy E^. Our analysis below shows that the first-order MIT is inevitably occur 
because of the nature of the long-range Coulomb interaction. Then we estimate the size of 
the transition for the case of fee YH3 based on the QSGVF calculation for insulator phase. 

Before looking into our model, let us remind the nature of the homogeneous electron gas 
Q, Q- At lower density n (density n can be specified by rg, 47rrg/3 = l/n), the Coulomb 
interaction v is more important than the kinetic energy. The exchange self-energy at E^ as 

^ = -1.2218/rs Ry dominates the kinetic energy e'^^'^ = |k^|V(2m) = 3.6832/r,2 

Ry at low density (large r^), where is the Fermi momentum. The Fermi energy is 
given as E^ = e"^™ + + S'^, where the Hartree term do not exist since it is cancelled by 
the background positive charge. These terms are evaluated at E^ and at for given n. 
Because of the behavior of e"^^"^ and S^, E^ as function of n is not monotonic. S^(n) oc —n^^^ 
dominates E^ for n — > 0, though e'^™ oc n^/^ does E^ for high n. S'^ enhances the effect 



of as seen in Table III in Ref. 
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by Hedin, where we see that E^ is negative at rg > 3. 
This large negative values of + S'^ overriding e"*™ at low n means the energy gain due 
to the Fermi statistics and the correlational motion of electrons (or the structure of the 
full many-body eigenfunctions). Since the Coulomb interaction is stronger at q ^ 0, the 
behavior of E^ is quite anomalous at n ^ 0. This is in contrast to models with short-range 
interaction, which is g-independent and we have oc —n. Then it can not dominate the 
kinetic term at low density. 

We treat a simplified model to avoid difficulty in real systems. In principle, we treat 
both of the insulator and the metal phases within QSGW{ot its extensions). However, in 
practice, simple sampling method for the Brillouin-zone summation [lO| is not applicable 
to the metal phase because we need to use too many k points to take into account the 
contributions from the small amount of holes in valence band (and electrons in conduction 
bands) . To avoid this problem, we consider a model with simplified energy bands with some 
assumptions. At first, we have to prepare for the insulator phase by QSGW. When we 
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can observe the MIT, should show small Eg, which is decreasing as a is being reduced. 
In our model, we only take into account the bands around the Fermi energy as shown in 
FiglU The valence bands are specified by the isotropic effective mass m^ai and its degeneracy 
A^vab the conduction band by rricon and A^^con- These parameters should be chosen so as to 
mimic the bands in the insulator phase given by QSGW. As we concentrate on the MIT, we 
only take into account the monotonic a-dependence in Eg] we assume mcon and mvai are not 
a-dependent. Further, we assume the rigid shift of energy bands (no deformations). 

For the model, we consider the situation that some amount of electrons (specified by 
density n) are moved from the valence bands to the conduction bands. Then we treat two 
quasi- Fermi energies, i?^^™'(n) for the electrons in valence bands, and £'^™'(n) for holes 
in conduction bands. As we discuss below, we will evaluate the total energy as function of 
n, and the energy minimum of the model occurs at some finite n (thus metal) below some 
critical Eg. 

-^vai™'(^) E^^^"^\n) are defined as the changing rate of the total energy per adding 
(subtracting) an electron. Thus the total energy AE{n) relative to the insulator phase is 
given by an adiabatic connection as 



where we can set E^^^'^'iO) = ^^f^r^O) = 0. Then is given as = e''"(^) + 



gHartrcc/^\ _j_ 5]x^^\ _|_ gg^^ ^^le elcctrou gas, e^™(n) = |k^P/(2mvai) oc 

The exchange part is given as = — 1.2218/(ers), where e means the effective dielectric 
constant representing the screening effect in the insulator phase, and is for the density 
n. We have to use n/N^ai instead of n if A^vai 7^ 1- We treat e as constant; e is little 
dependent on a in the case of indirect gap as in fee YH3. For S^, we use the RPA formula, 
e.g, see Eq.(86) in Ref . jisl . where we use W = v/e{q, uj) = v/e x 1/(1 — vxo/e)- This W is 
obtained in the RPA with v/e instead of the bare Coulomb interaction v. xo here contains 
the contribution due to electrons in valence bands and due to holes in conduction bands, 
but no interband contributions; neglecting the interband contribution will be reasonable for 
the case with indirect gap. We neglect g^'^''*''*^*^^ because g^^''*'''^*^ oc n can be neglected in 
comparison with at least for small n. We further neglect the other kind of correlational 
effect beyond QSGVFas excitonic effects between electrons. With these assumptions, we can 
evaluate i?^™'(n) and also E^^l™\n) for given mvai, A^vai, ""^con, A^con, and e. In our model, 
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a-dependence is only in Eg] the dependence cause just a constant shift in the integrand of 
Eq. ([1]) as function of n. 
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FIG. 2: Solid line shows — £'^^™(n) for the case rrivai = "icon = unit of electron 

mass), A^vai = -^con = 1) and e = 8. Dotted line is without T,'^ contribution. The difference of areas 
A— B shows the total energy difference between insulator and metal. When Eg < 0.08 eV, we have 
B > A resulting the stability of the metal phase at Q. See text. 



Let us look into the numerical results. In Fig 121 solid line shows 
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1 (in unit of electron mass), A^vai = 



of the electron-hole symmetry in this case, E^^^ 
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1, and e = 8. Because 
2EEZ'^Hn) is satisfied. 



Derivative of the line is divergent at n — since it is oc —n^^^. At high density, its behavior 
is controlled by e'^'" oc n^^^. The integrand Eg + i?^™'(ra) — E^^i'^^n) in Eq. ([1]) is shown 



by the solid line {—Eg as zero level). We show the case with Eg = O.OSeV. At the crossings 



P and Q, we have Eg + E'^Z'^^n) - E^^, 



Fermi , 



n) 

(fAE 



0, that is, the derivative 



0. Q 



corresponds to the stable solution satisfying > 0. The energy AE of the metal phase 
at Q is given as the difference of areas A— B. We have A— B=0 at Eg = O.OSeV in the case; 
for E^ < El 



.g O.OSeV, the metal phase at Q has lower energy than the insulator phase 
since A<B. At Q point the band overlap is E'^^^ =0.09 eV, where we simply assume the 
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TABLE I: Calculated size of MIT in our model, for given e, mvai, mcon-, -^com (we use A^vai = !)• 
When Eg is getting smaller and reach at Eg = E"^ , it cause the first-order phase transition to 
metal with the band overlap E^^^- 
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rigid band shift (no band narrowing nor widening). The dotted hne is without S^; we see 
that the contribution from Tf^ enhance E'^ almost twice larger. Considering the behavior of 
-^co™(^) ~ -^vaf'^H'^)' this first-order phase transition can be a general phenomena for the 
band-overlap MIT. 

In Table [U calculated and E'^^-^ are given for kinds ^vai, '^con, -?Vvai and e. As mvai 
and mcon are heavier, E'^ and E'^^^ g^t larger because the contribution from kinetic term 
is smaller, e, which determines the size of the effective Coulomb interaction, can strongly 
affect on E'^ and i^^^j. We may overestimate S'^ a little in RPA, thus we does E'^; this is 
indicated by the fact that for homogeneous electron gas given by the accurate method 
[isl is 69 % of RPA value at rs=4 (62 % at =100). The parameter set m^^\ = .41, mcon = 
.92, Neon = 4, and e = 8 is to mimic the energy bands of fee YH3. Its bottom of conduction 
band is at L point, where two Fermi surface exist, but we neglect the smaller one because it 
contains only ~10% of electrons of bigger Fermi surface. The Fermi surface is anisotropic, 
but we take the simple average as mcon = {fn^mymr^Y^^ = 0.92. The valence top is at F, 
where it has one spherical Fermi surface. We will report how the band gap changes as the 
function of a elsewhere together with other analysis for comparison with experiments [l^. 
In anyway, our result of E^'^ = O.OSeV should be taken as a semi- quantitative prediction 
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since out model treatment is very simplified. 

Though our MIT mechanism may occur in reality, we have not found experiments which 
are directly related to our MIT mechanism. As we discuss elsewhere 16|] , the MIT for 
YH3 observed in experiments will be mainly controlled by the structural transition, thus 
the band-overlap mechanism here will be not directly related to the experiments. In the 
case of GdN [l7|, the MIT explained here might occur. See Fig.2 in RefilTI. However, we 



need to improve our treatment to have some numerical prediction for GdN because it is 
magnetic and multiple bands are involved. As for the bcc hydrogen [sl, we may apply our 
MIT mechanism. In fact, the variational Monte-Carlo results seems to indicate the weak 
first-order phase transition 18|, [lOj . 

In conclusion, we have given a theoretical analysis for the metal-insulator transition 
through the band-overlap mechanism. We have showed that it occurs as the first-order 
phase transition; when the band gap is getting smaller, the insulator phase suddenly changes 
to the metal phase because of the energy gain due to the Fermi statistics and due to the 
correlational motion for the Coulomb interaction. In our treatment, the self-consistency in 
QSGVFis not fully included; we only consider the rigid shift of energy bands. If we fully 
include the self-consistency, we should have some deformation of the energy bands (changes 
of effective masses). However, we expect that our conclusion will be unchanges qualitatively. 
To treat kinds of materials in practice, it will be necessary to improve our method. Further, 
we will have to examine roles of other effects like excitonic effects or phonon effects. 

This work was supported by ONR contract N00014- 7-1-0479. We are also indebted to 
the Ira A. Fulton High Performance Computing Initiative. 
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